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ABSTRACT 



v~j ■ Massive tori with w 1O 8_9 M are predicted to extend on ~100 pc scale around 

the centre of elliptical galaxy progenitors by a model of a supermassive black hole 
(SMBH) growth coeval to the spheroidal population of the host galaxy. Direct 
detection of such massive tori would cast light on a key physical condition that 
allows the rapid growth of SMBHs and the appearance of QSOs at high redshift. 
For this reason, we examine the detectability of such structures at substantial 
redshift with the Atacama Large Millimeter Array (ALMA). We propose that 
submillimeter galaxies (SMGs) are the best targets to test our predictions. In 
order to assess the observational feasibility, we estimate the expected number 
counts of SMGs with massive tori and check the detectability with the ALMA 
instrument, the unique facility which can resolve the central region of high red- 
shift objects. Our work shows that ALMA will be able to resolve and detect 
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high-J (J > 4) CO emissions from ~100 pc scale extended massive tori up to 
z ~ 2. Observations of lensed SMGs will yield excellent spatial resolution, al- 
lowing even to resolve their massive tori at higher redshift. We discuss further 
the detectability of the HCN molecule, as a better tracer of the high density gas 
expected in such tori. The final goal of these kind of observations is to pinpoint 
possible physical mechanisms that storage in the very central galactic regions 
very large amount gas on timescale of several 10 8 yr. 

Subject headings: galaxies:nuclei — galaxies:evolution — galaxies:formation 

1. Introduction 

The energy emitted by active galactic nuclei (AGNs) is commonly ascribed to accretion 
onto a supermassive black hole (SMBH). Recent high- resolution observations of galactic 
centers indicate the presence of a SMBH whose mass correlates with the mass, the velocity 
dispersion, and the luminosity of the galactic bulge (e.g., Kormendy & Richstone 1995; 
Richstone et al. 1998; Magorrian et al. 1998; Laor 1998; Ferrarese & Merritt 2000; Tremaine 
et al. 2002; McLure & Dunlop 2001; 2002; Marconi & Hunt 2003). QSO hosts are mostly 
luminous and well evolved early-type galaxies (e.g., McLeod & Rieke 1995a; Bahcall et 
al. 1997; McLure, Dunlop & Kukula 2000; Dunlop et al. 2003). Taken together, these 
observations indicates that the formation of a SMBH, a bulge and a QSO is closely related 
to each other. While it is still not clear which is the main mechanism for the formation 
of a SMBH, many studies suggest that a BH growth by gas accretion processes plays an 
important if not dominant role (Salucci et al. 1999; Fabian & Iwasawa 1999; Yu & Tremaine 
2002; Hosokawa 2002; Marconi et al. 2004; Shankar et al. 2004). 

Kawakatu, Umemura & Mori 2003 (hereafter KUM03) and Granato et al. 2004 (here- 
after G04) proposed a physical model for the co-evolution of QSO BHs and the spheroidal 
component of galaxies. The formation of these latter is computed with a semi-analytical 
technique, which includes the effects of stellar and QSO feedback, and model outcomes have 
been more widely compared to observations of several galaxy populations (G04; Silva et al. 
2005; Granato et al. 2006; Lapi et al. 2006). In these models the fueling of the central BH 
in forming spheroids takes place in two steps: (i) a low-angular momentum gas reservoir is 
formed at a rate essentially proportional to the star formation rate on a time scale of a few 
10 8 yr, then (ii) the reservoir gas is accreted onto the central BH at a rate of the order of 
the Eddington limit. Since in the earlier phases of a SMBH growth the former rate largely 
exceeds the latter, a massive reservoir accumulates; most of its mass is accreted onto the 
SMBH only at the last e-folding times (t e ~ 4 x 10 7 yr) and powers an enormous energy 
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release, which stops both the star formation and the growth of the reservoir. The massive 
reservoir is likely to be a geometrically thick massive disk (we call it massive tori) because 
of some residual angular momentum. We will comment on the passage from the massive 
reservoir to the massive torus in details later. A plausible physical mechanism responsible 
for step (i) is the radiation drag (Umemura 2001; Kawakatu & Umemura 2002), thus here- 
after we refer them as the radiation drag model. It should be however kept in mind that any 
other mechanism producing a similar rate of production of massive tori, would not alter the 
essence of these computations. 

The scenario proposed by KUM03 and G04 establishes a well defined sequence con- 
necting various populations of massive spheroidal galaxies: (i) vigorous and rapidly dust- 
enshrouded star formation activity, during which a central SMBH grows, and explaining 
the sub- mm galaxy population (SMG), which are examples of high-z ultraluminous infrared 
galaxies (ULIRGs); (ii) QSO phase halting subsequent star formation and (iii) essentially 
passive evolution of stellar populations, going through an Extremely Red Object (ERO) 
phase (see Fig. 1; For details Fig. 2 and Fig. 3 in KUM03 and Fig. 3 in G04). 

Here we suggest that the massive torus, before accreting onto the SMBH (time-scale is 
« 10 8 ~ 9 yr), forms a massive torus of the kind envisaged by unified models of AGNs (e.g., 
Antonucci 1993; Granato & Danese 1994; Urry & Padovani 1995). Thus the radiation drag 
model predicts the existence of a massive torus with « 1O 8_9 M in the early phase of SMBH 
and spheroid growth (stage (i) and also the beginning of stage (ii), see Fig. 1). The impact 
of the direct detection of massive tori is to reveal whether the large mass ratio of a torus to 
a black hole, M t0 rus/^BH ^> 1, is a key physical condition which triggers the rapid growth 
of SMBHs and the appearance of QSOs. This issue will be discussed in details later. Hence, 
we investigate the required spatial resolution and expected CO and HCN emission lines from 
massive tori. We focus our investigation on the high- z universe (z > 1) where the peak 
of QSO activity occurs (around z ~ 2) and mild AGN activity has already been detected 
in high- ,2 SMGs (Alexander et al. 2003, 2005a, & 2005b) in very good agreement with the 
G04 predictions (Granato et al. 2006). In particular, we here investigate the detectability of 
massive tori at z > 1 with the next generation millimeter interferometer, ALMA (Atacama 
Large Millimeter Array). The rest of this paper is organized as follows. In §2, we briefly 
review the formation of massive tori in the early phase of the BH growth and set up the 
model for massive tori. In §3, we first estimate the expected spatial resolution and expected 
CO and HCN fluxes from massive tori. Secondly, we discuss the best targets with massive 
tori and then show the number count of candidates with massive tori. Finally, we mention 
the importance of the velocity mapping on massive tori, and the impact on the detection 
of them for SMBH growth. Our discussions and conclusions are given by §4. Through this 
paper, we adopt the cosmology indicated by WMAP data (Bennett et al. 2003; Spergel et 
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al. 2006), i.e., a flat Universe with the Hubble parameter H = 70kms 1 Mpc 1 , f2 M = 0.27 
and tt A = 0.73. 



2. Model 

2.1. Mass accretion due to radiation drag 

We review a physical model of the formation of a SMBH in forming spheroids, based on 
the radiation drag-driven mass accretion. The radiation drag, driving the mass accretion, 
originates in the relativistic effect in absorption and subsequent re-emission of the radiation. 
This effect is naturally involved in relativistic radiation hydrodynamic equations (Fukue, 
Umemura, & Mineshige 1997). The angular momentum transfer in radiation hydrodynamics 
is given by the azimuthal equation of motion in cylindrical coordinates, 

where E is the radiation energy density, is the radiation flux, P^ is the radiation stress 

tensor, and Xd is the mass extinction coefficient which is given by Xd — n&o&j p gas with 

the number density of dust grains n^, the dust cross-section <7& and the gas density p gas . By 

solving the radiative transfer including dust opacity, Kawakatu & Umemura (2002) evaluated 

the radiative quantities in a clumpy medium, E, F^, and P^, and thereby obtained the total 

angular momentum loss rate. Then, these authors estimated the total mass of the dusty ISM 

accreted on to a central massive reservoir, M rcs . By using the relation, M gas /M gas = — J/J, 

where J and M gas are the total angular momentum and gas of ISM. Then, the res ervoir mass 

is assessed as t t ' 

M TCS (t) = [ M gas dt = - [ M gas -dt. (2) 
Jo Jo J 

In the optically-thick regime of the radiation drag, — M gas J/J = L sp h(t)/c 2 , where L sp h(t) 
is the total luminosity of spheroids. The radiation drag efficiency depends on the optical 
depth r in proportion to (1 — e~ T ) (Umemura 2001). Thus, the total mass of reservoir can 
be expressed by 

M ms (t) = % rag jf (l - e-^')) dt>, (3) 

where r(t) is the optical depth of spheroids measured from the center. Here, we estimate 
the evolution of r{t) by using an evolutionary spectral synthesis code 'PEGASE' (Fioc & 
Rocca-Volmerange 1997). The efficiency r/drag is found to be maximally 0.34 (Kawakatu & 
Umemura 2002). 
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Moreover, it is worth keeping in mind the ratio of reservoir mass to spheroidal mass 
(see Umemura 2001 in details). The radiation energy emitted by a main sequence star is 
0.14e times the rest mass energy of the star m^c 2 , where e is the energy conversion efficiency 
of nuclear fusion from hydrogen to helium, which is 0.007. Thus, the total luminosity of 
spheroids is estimated to be 

L sph (t) ~ 0.14 e 5 sph (^c 2 , (4) 

where S , sp h(t) is the star formation rate in the spheroidal system. By substituting this in 
equation (3), the reservoir mass to spheroidal mass is given by 

M res(t) _. n w 1A - 3 



M sph {t) 



2 X 10~ d , (5) 



which is comparable to the observed black hole-to-bulge mass ratio in nearby galaxies. Here 
M sph = f*S sph (t')dt'. 



2.2. Formation of a massive torus 

In §2.1, we introduced the res ervoir mass in the context of the radiation drag-induced 
mass accretion. However, the massive reservoir does not evolve into a SMBH directly, 
because the radiation drag can not remove the angular momentum thoroughly (Sato et al. 
2004). Thus, some residual angular momentum would terminate the radial contraction. 
In the massive reservoir, the viscosity is expected to work effectively because the viscous 
timescale is shrinked by the radiation drag (Mineshige et al. 1998). Thus, the massive 
reservoir is likely to be a massive self-gravitating viscous disk. In such a massive disk the 
extreme star formation occurs in the outer part of a massive torus via the gravitational 
instability, since Toomre Q parameter is less than unity, where /tc s /7rGE g for the epicycle 
frequency k, the sound velocity c s and the surface density S g . If this is the case, the massive 
res ervoir will become a geometrically thick massive disk (hereafter we call it a massive torus: 
M torus = M res ), by supporting the vertical thickness by the energy feedback from supernovae 
(Wada & Norman 2002) and radiation pressure from the starburst (Ohsuga & Umemura 
2001; Thompson et al. 2005; Watabe & Umemura 2005). 

As for the BH growth, we here assume that the mass accretion driven by the viscosity 
onto the BH horizon is determined by an order of Eddington rate, Msdd = L Edd /c 2 ~ 
O.2M0yr _1 (MBH/lO 8 M ), where L E dd is the Eddington luminosity. As mentioned in §2.1, 
the mass accretion rate from a host galaxy to a galactic center via radiation drag is Md rag ~ 
L sph /c 2 ~ O.l(L sph /lO 12 L )M yr- 1 . According to KUM03 (also G04), the starburst phase 
is closely related to the formation of a QSO. During the first stage, vigorous star formation 
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occurs, while a SMBH is still growing (M BH ~ 10 6 — 1O 7 M ). In this phase, the existence of 
the massive torus with M torus /M B H 3> 1 can be predicted, since Md rag largely exceeds M Edd . 
At the end of starburst phase, the SMBH has grown and radiates at maximum luminosity 
(QSO phase). In the QSO phase, M BB approaches M torus and thus M torus /MBH — 1- After 
the QSO phase, the AGN luminosity decreases rapidly, since almost all of the matter in 
massive tori has fallen onto the central BH, namely QSO would evolve into low luminosity 
AGNs (LLAGNs). Hence, M torus /M BH can be expected to be much smaller than unity 
for spheroidal galaxies with LLAGNs. This is consistent with observation which elliptical 
galaxies and galactic bulges have no gas. Therefore, it can be concluded that the evolution 
of the mass ratio (M torus /M B H) reflects the history of the SMBH growth. 



2.3. Physical parameters of massive tori 

The radiation drag mechanism, as mentioned §§ 2.1 and 2.2, predicts the presence of 
massive tori endowed with a mass M torus ~ 10 8 — 10 9 M equal to the final mass of BHs 
powering QSOs (see Fig.l). Typical length scales of such tori can be expressed by the 
accretion radius of the "BH plus massive torus" system (r a ). The accretion radius is given 
by r a = C(M B h + -^torus)/ cr sph> where a sp h is the velocity dispersion of spheroids. Following 
the virial theorem for spheroids (of h = GM sph /r h ), the radius can be written as 

M BH + M torus ^ /M torus \ 

r b - -T7 — r b; (6) 



M sph J \ M sph j 

where r b is the effective radius of spheroids. Here we consider the case of M BH <C M torus 
since we focus on the rapid growth phase (see stage (i) in Fig. 1). By substituting equation 
(5) in equation (6), the r a can be estimated by r a ~ 20 pc (r b /10kpc). In the present 
paper, the size of a massive torus is defined as r tor us = Ar a . The factor A can be inferred 
by high resolution observations of tori in nearby AGN (e.g., Jaffe et al. 1993; van der Marel 
et al. 1998; Davis et al. 2006) and turns out to be at least A R3 5. Recent IR and X-ray 
observations support ~100 pc scale extended massive tori (e.g., McLeod & Rieke 1995b; 
Maiolino et al. 1995; Maiolino & Rieke 1995; Granato et al. 1997; Malkan, Gorjian & Tran 
1998) rather than the compact massive tori. However, r t0 rus (or A) is physically determined 
by the quantity of residual angular momentum removed by the radiation drag process. In 
order to reveal this issue, we should explore the sophisticated numerical simulations, which 
are our future study. Thus, a compact massive torus can not be ruled out physically. In this 
paper, we examine two cases; one is a compact torus with r toms = 20 pc (A=l), the other is 
an extended torus with r t0 rus = 100 pc (A— 5). 



Concerning the rotation law and the surface density profile in massive tori, we assume 



-7- 



v<t>(r) = v^ (r/r torus )- a and S g (r) = S gi o(r/r t orus)" /3 , where v^o and S gj0 are the quantities 
at r toms . S gj0 is related to the torus mass through M torus = J^ tOTUS 27rrS g (r)dr, where r in is the 
sublimation radius for graphite grains. The size r in is expressed as r in = 1.3 L^y 46^1500 P c > 
where Luv,46 is the AGN UV luminosity in unit of 10 46 erg s -1 , and T 150 o is the grain sublima- 
tion temperature in unit of 1500 K (Barvainis 1987). Recently, Suganuma (2006) measured 
the inner edge of dusty torus r- m by the reverberation mapping technique, and thier ob- 
servational results strongly support that r in is determined by the dust sublimation radius. 
Although the size r in depends on the size of dust grains and dust composition (e.g., Laor & 
Draine 1993), these does not affect on the expected line flux significantly as far as r in <C r tor us 
. Thus, we set up r in =l pc in the present paper. The rotation velocity, v^o = A~ a v l j ) {r- 3 ), 
can be determined from t></,(r a ) = <7 sp h. Since a and (3 are unknown free parameters, we will 
calculate the CO flux density corresponding to the following parameter set (a = —1,0,0.5 
and = 0,1,2). 

Provided that the scale height (/itorus) of massive tori is ~0.3xr tO rus (e.g., Wada & 
Norman 2002), the average number density and column density of molecular hydrogen are 
ra(H 2 ) « lO 3 cm- 3 (M gas /lO 8 M )(r torus /lOOpc)- 3 andiV(H 2 ) « lO 23 cm- 2 (M gas /lO 8 M )(r torus /lOO pc)" 2 , 
respectively. Thus massive tori are optically thick and thermalized at CO lines. This happens 
to be for any molecular clouds with hydrogen column density larger than iV(H2) > 10 21 cm -2 
(e.g., Tielens & Hollenbach 1985). For /i tor us, we expect that the opening angle is almost zero 
in the quite extreme conditions we are envisaging. However, the expected CO flux density 
does not depend on h torus as long as the massive tori are optically thick for CO lines. 

As for the gas temperature T g of the optically thick massive torus, Ohsuga & Umemura 
(2001) analyzed it by using one-dimensional radiation hydrodynamic equations, coupling 
with ionization processes and thermal processes. They found that dust cooling is effective 
and thus the gas temperature T g is w 100 K. Thus, we assume T g = 100 K in this paper. This 
is also supported by the recent numerical results that the massive tori would be dominated by 
the cold (T g ~ 100 K) and dense (n(H 2 ) ~ 10 2 cm~ 3 ) molecular gas (e.g., Wada & Tomisaka 
2005). Finally, we note that the CO luminosity and the CO-to-H2 conversion factor are a 
function of the metallicity (Radford, Solomon & Downes 1991) and the far-UV flux from 
AGNs and/or massive stars lower of a few magnitudes for n(H 2 ) > 10 2 cm~ 3 (e.g., Tielens 
& Hollenbach 1985; Mochizuki & Nakagawa 2000). 
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3. Results 

3.1. Required spatial resolution 

Fig. 2 shows the spatial resolution required to detect 100 pc scale extended massive 
tori. A resolution of ~ 0.01" is necessary, well beyond the capabilities of present-day mm and 
sub-mm interferometers (typically ~ 0.5"). By converse, the angular resolution of ALMA 
from 86 - 720 GHz with 9(") ~ 0.005(")(z/ obs /600 GHz)" 1 0, is sufficient to resolve high- 
J (J > 4) molecular emissions from 100 pc scale extended massive tori, where is the 
observed frequency. On the other hands, 20 pc scale compact massive tori cannot be resolved 
unless the objects are gravitationally lensed (e.g., Downes and Solomon 2003). The lensing 
effect magnifies not only the intrinsic CO line luminosity, but also the scale of emitting 
region. If we assume the magnified CO image by an ellipse, the magnification factor is 
/i = a/r torns , where a = .DA©obs/2 is the apparent semi-major axies of the magnified image, 
where Da and o b s are the angular distance and the observed angular length of magnified 
image, respectively. Thus, it is possible to achieve the super-resolution which can be \x times 
as high as the instrumental resolution, e.g., \i = 2 — 50 for lensed SMGs (Solomon & Vanden 
Bout 2005). As seen in Fig. 2, for lensed objects we can resolve lower J-CO lines from 
massive tori and/or study massive tori at higher redshift. 

3.2. Expected flux of the CO molecule 

Since most of molecular clouds in massive tori are in cold and dense (§2), it is expected 
that CO molecular lines are good tracers of molecular gas mass. In this paper, we express 
the CO line luminosity l' co (J, J — 1) in unit of Kkms" 1 as the product of the velocity- 
integrated source brightness temperature T b Aw rest , where Au rest is the rest-frame line width. 
Then, the CO line luminosity in the arbitrary transitional level from J to J — 1 is given 
by l'co(J, J — 1) = TbAuregt- Thus, the CO line luminosity ratio for two lines, Rj t j-\ = 
l' C0 {J, J — 1)//q O (1,0), in the same source is equal to the ratio of the intrinsic brightness 
temperature. Since we deal with the case that CO lines are optically thick and thermalized 
up to J = 6, we set Rj,j-i = 1. Note that the gas temperature can be warmer than in 
normal galaxies due both to the intense radiation field via the huge star-forming activity 
and to the accretion onto their central growing BH. A flat distribution of Rayleigh- Jeans 
brightness temperature up to J=6 has been reported in the starburst galaxies (e.g., Kawabe 



1 This comes from the ALMA webpage (http://www.eso.org/projects/alma/science/bin/sensitivity.html) 
See note 3 in this webpage. 
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et al. 1999). 

The line flux density /co(<A J ~ 1) i n um t °f Jypc -2 corresponding to the CO line lumi- 
nosity l' C0 (J, J — 1) can be derived. Using H 2 mass-to-CO line luminosity relation l' co (J, J — 
l)(r) = i?j j j_iS(H 2 )(r)/Xco, the corresponding received CO flux density fco(J, J — 1) is 

2k ^ (J, J - l)(r)z/ 2 cst (J, J - 1)(1 + g) cosji) 
fco(J,J-l)(r)-- ^—^ w , (7) 

where k is the Boltzman constant, % is inclination angle of the massive tori, the fco{J, J — 1) 
is the observed flux density in Jypc -2 , S(H 2 )(r) = S g (r), At> rest (r) = tfy(r), ^ re st(^ </ _ 1) = 
115 J GHz is the rest frequency of the transition and is the luminosity distance. Here 
Aco = O.8M (Kkms" 1 )" 1 is a typical value adopted for ULIRGs and SMGs (Downes 
& Solomon 1998; Solomon & Vanden Bout 2005). This value is comparable to Xco — 
0.5M Q (Kkms -1 ) -1 , which was calculated by three-dimensional, non-local thermal equilib- 
rium radiative transfer calculations for high- J CO lines (e.g., J = 3 and J = 4) emitted 
from inhomogeneous dusty tori (Wada & Tomisaka 2005). 

The integrated CO flux density (Jy) of massive tori over their surface can be calculated 
by S co (J, J — I) — 27rr/ co (J, J - l)(r)dr. Then, 

Sco(J,J-l) « 10 3 Jj J 2 H(a = 0.5, P = 1) (^y 1 (^f^ 

M torus \ / Av rcstfi \~ 1 + z 



3 x lO 8 M y V100km/s J D 



where Av Testfi = lOOkm/s (A/5)~ (a/a5) (M torus /3 x 10 8 M o ) a5 (r a /20 pc) -0 5 , a inclination an- 
gle being equal to 45 degree, and H(a,/3) is the function of a and (3. We summarize the 
value of H(a,/3) in Table 1. 

Fig. 3 shows the expected CO emission from 100 pc scale extended massive tori as a 
function of redshift for a fiducial case (M torus = 3 x 1O 8 M , A = 5, a = 0.5, and (3 = 1). 
The five lines denote the CO flux density Sco(J, J — 1) at the transition level (J, J — 1). 
At the left side of filled circles, each line can be resolved with ALMA (see Fig.l). The 
CO lines without filled circles will not be resolved. The horizontal line shows the ALMA 
detection limit (50 antennas) 0.5 mJy requested 5a in 25 km/s channel in 12 hours of on- 
source integration time[§. The high- J CO lines (J > 4) can be detectable and resolvable up 
to z = 2 with ALMA. Again it is worth recalling that for gravitational lensed objects we can 



2 Note that the ALMA detection limit for lines depends on the observed bands as follows; 
0.31 mJy for band 3 (86-116 GHz), 0.42 mJy for band 4 (125-163 GHz), 0.39 mJy for band 5 
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detect lower J-CO lines and/or we can study higher redshift objects, because the observed 
CO flux density increases /i times with respect to that estimated in eq.(8). 

Finally, we discuss the dependence of A (or r torus ), a and /3 for Sco(J, J — !)■ First, we 
investigate the effect of the function H(a,(3). As seen in Table 1, it is found that H(a,/3) 
is order of unity, except for the case of a = —1 and [3 = 2. In the case of a = —1 and 
[3 = 2, the expected CO flux density Soo(J, J — 1) has larger value compared with other 
cases, because of a high gas concentration and slower rotation velocity in the inner region 
of massive tori. As mentioned in §2.3, Af res t,o is linked to the size of the torus. Thus, 
the emissivity Sco(J, J — 1) also depends on the difference of r torus (or A). Then, the CO 
flux density from compact massive tori (A = 1, a = 0.5 and < [3 < 2) is ~ 3 times 
as small as fiducial case, because Af rest0 of compact massive tori is larger than that in 
the fiducial case. Thus, the detection of compact massive tori with ALMA would be more 
difficult than the extended massive tori, taking into account of the spatial resolution (§3.1). 
Thirdly, the variation of a also affect on Sco{J> J ~ 1) f° r the extended massive tori (A = 5 
and < (3 < 2). From eq.(8), Sco(J, J — 1) of extended massive tori with rigid rotation 
(a = —1) and flat rotation (a = 0) are about 0.1 and 0.5 times smaller that in the fiducial 
case, respectively. From this, we find that it will be difficult to detect the extended massive 
tori with ALMA if they rotate rigidly. To sum up in §3.1 and 3.2, we conclude that 100 pc 
extended massive tori can be resolved and detected with ALMA, except for the extended 
rigid rotating massive tori. 

3.3. Expected flux of the HCN molecule 

HCN emission was observed locally in giant molecular clouds and in nearby galaxies 
(see i.e. Gao and Salomon 2004 and references therein). At high redshifts only gravita- 
tionally lensed sources have been detected, namely H1413+117 (Cloverleaf) at z=2.Q, IRAS 
F10214+4724 at z=2.3, and the quasar APM08279+5255 at 2=3.911 (Solomon et al, 2003; 
Vanden Bout et al., 2004; Wagg et al, 2005). 

To estimate the flux of the high- z targets discussed here, we assume that the HCN 
molecule has the same velocity width of the CO line, since both emissions are originated 
within the same dynamical structure (Wagg et al., 2005). As observed for CO lines, we 



(163-211 GHz), 0.46 mJy for band 6 (211-275 GHz) and 0.50 mJy for band 7 (275-370 GHz), 
wchih can be obtained by the ALMA sensitivity calculator as mentioned on the ALMA webpage 
( http: / / www.eso.org/ projects / alma/science /bin /sensitivity.htm| . Here we take the greatest value (0.50 
mJy). 
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suppose also for HCN luminosity a proportionality with the molecular hydrogen mass. For 
normal galaxies a relation between CO and HCN luminosities is established (Gao & Solomon, 
2004), typically ZhcnAco ~ 0.1. For most ULIRGs with AGNs this ratio turns to be higher, 
> 0.2 (Imanishi et al. 2006). We take here the value of 0.2, similar to that of Mrk 231 
and IRAS 17208-0014, and not far from that of APM08279+5255 (0.34). The HCN line 
flux is then /h CN (1,0) = 0.3(XHCN/3.3) _1 (Xco/0.8) _1 S g , which is similar to the relation 
^hcn(1>0) = Sg/XncNj where X H cn ~ 7 M (Kkms _1 ) for warm gas (Wagg et al. 2005). 

In Fig. 4 we compare the expected HCN line emission from 100 pc scale extended 
massive tori (Mbh = 3 x 10 8 M q , A = 5, a = 0.5 and /3 — 1) as a function of redshift to 
the ALMA sensitivity limit. The dependences of A (or r toms ), a and (3 are same as the case 
of CO lines (§3.2). The detection of HCN lines from massive tori is more difficult because 
of the lower abundance of this molecule with respect to CO. At z « 1 the detection of the 
central massive torus in this line would require more than 24 hours of on-source integration 
in the highest angular resolution configuration. In lower resolution configurations this line is 
detected at any redshift, however, the constraint on the model presented in this model will 
not be stringent. It will be possible to detect and resolve HCN lines for gravitational lensed 
objects. 



3.4. Which objects are associated to massive tori ? 

We discuss here the best targets associated with massive tori and their expected number 
counts. Observationally, the majority of SMGs is coeval to QSO activities with z ~ 2. From 
their properties (the massive dynamical mass, the high star formation rate, etc.), SMGs 
are considered to be progenitors of present-day elliptical galaxies (e.g., Smail et al. 2004). 
Recent ultra-deep X-ray observations have suggested that > 50% of SMGs have AGNs whose 
luminosities are relatively lower than those of coeval QSOs (Alexander et al. 2005a). Thus, 
following our reference model we expect that massive tori are already in place in these 
galaxies and they are feeding the growth of the central BH with M B h ~ 10 6 — 1O 7 M (Borys 
et al. 2005). Recently, Menedez-Delmestre et al. (2007) and Valiante et al. (2007) observed 
the mid-infrared spectra of a sample of SMGs, in which they found that in the rest frame ~ 
6-8 /iiii there is a contribution from a relatively weak AGN continuum emission to be added 
to the contribution from star formation (cfr Fig.3 in Valiante et al.'s paper). Thus, their 
results stress once again that the most interesting targets in the search for tori could be 
the submm galaxies with 24 /xm counterparts. According to KUM03 (also G04), the dust 
emission from massive tori in SMGs is cooler than that in QSOs because of the higher column 
density of massive tori and relatively weaker AGN activity. It will be also worth comparing 
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the rest frame infrared color (e.g., f(25 /xm)/f(60 /iin)) of SMGs with that of coeval QSOs. 
Moreover, Greve et al. (2005), Tacconi et al. (2006) and Gao et al. (2007) have reported 
~10 CO detected SMGs and ~ 5 HCN detected SMGs, indicating the presence of large 
masses of dense gas (~ 1O 1O M ) within ~ 4kpc. SMGs are largely present among the Early 
Universe Molecular Emission Lines Galaxies (EMGs) listed by Solomon and Vanden Bout 
(2005); a significant fraction are also lensed objects. 

Therefore, SMGs (especially SMGs with AGNs) are good candidates to explore a massive 
torus and to test the radiation drag model. We evaluate the expected number count of 
SMGs having massive tori with the G04 model. Note that the G04 model predicts number 
counts, redshift distribution, and AGN accretion rates of SMGs in very good agreement with 
observations (G04, Silva et al. 2005, and Granato et al. 2006). Fig. 5 shows the redshift 
distributions of bright SMGs with different masses of massive tori. The expected number 
count of galaxies with M torus > 10 8 M Q is ~ 100 — 200 deg~ 2 at 1 < z < 3. By comparing our 
predictions with ALMA data, it would be possible test and refine the radiation drag model. 

Lastly, scaled down versions of star forming galaxies harboring an AGN can be studied 
in the local Universe. The narrow line Seyfert galaxies (NLSls) and ULIRG (especially 
type I ULIRGs) are characterized by a smaller BH and a higher mass accretion rate (e.g., 
Pounds et al. 1995; Boiler et al. 1996; Mineshige et al. 2000; Mathur et al. 2000; Collin 
& Kawaguchi 2004, Hao et al. 2005; Kawakatu et al. 2006; Kawakatu et al. 2007). Thus, 
they would correspond to the early phase of the BH growth. If massive tori are detected 
in targets at z < 0.01, they will be easily resolved with present mm interferometers as 
well. Indeed, CO emissions from < 100 pc has already detected in nearby NLS1 NGC 4051 
(L. Tacconi, private communications). By investigating the existences of massive tori in 
SMGs (spheroids), ULIRGs (spheroids and disk galaxies) and NLSls (disk galaxies), it will 
be possible to clarify whether the BH growth process depends on the morphology of host 
galaxies and redshift. 



3.5. Velocity mapping on massive tori 

For low- ,2 targets, ALMA can study not only the existence of massive tori, but also the 
sub-pc structures of massive tori, although it is hard to do it for high-z objects. Thus, we 
will make the channel maps of molecular lines with ALMA, as Davis et al. (2006) carried 
on it for nearby Seyfert galaxies NGC 3227. As a consequence, it is possible to distinguish 
between the circular velocity, and the velocity dispersion, a in the region above the spatial 
resolution of ALMA, that is, from a few pc to ~ 100 pc (e.g., Fig.3 in Wada and Tomisaka 
2005). If a is comparable to in massive tori, it means that the kinetic energy is dominated 
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by the random motion. In order to maintain the vertical structure the continuous energy 
injection should be necessary because the cooling time of molecular clouds via collisions is 
comparable to the orbital period. As mentioned in §2.2, the radiation pressure via starburst 
and/or the internal turbulence caused by supernova explosions may be a plausible origin of 
its energy source. If the scale height of the torus is determined by a balance between the 
vertical components of the centrifugal force due to a central SMBH and the vertical velocity 
dispersion (e.g., Wada & Norman 2002), the ratio a/v^ reflects the ratio /itorus/^torus; where 
^toms is the scale height in the body of the massive torus. Hence, the velocity mapping of 
massive tori is helpful in providing a quantitative estimate of the origin of geometrically 
thick obscuring torus. 

In addition, it is expected that the starburst luminosity in a massive torus is closely 
related to its vertical structure and the determination of this through ALMA observations is 
a fundamental test of this model. ALMA velocity mapping assesses the velocity dispersion 
a and the scale height /i tor us for low-z massive tori and thus it allows to estimate the viscous 
coefficient of turbulent motion z/ tu b ~ o"/itorus, assuming the largest eddy size is comparable 
to fotorus- Using z/ tu t>, we can estimate the mass accretion rate due to the turbulent viscosity 
M B h ~ ^torus/^vis where t vis ~ r 2 /u tuh . Thus, it will be able to check whether the mass 
accretion rate via the turbulent viscosity is a key element in explaining the AGN activity, 
because both M torus and i/ tu b are observable. In order to get an insight into the physics 
of a fueling process, it would be valuable to derive the radial velocity v T (r), which would 
depend on the mass accretion process (e.g., the gravitational torque from a non-axisymmetric 
potential, the turbulent viscosity, the radiation avalanche, etc ). Therefore, the velocity 
mapping on massive tori is essential to explore the formation of obscuring massive tori and 
also the fueling process from ~100 pc to ~1 pc, which is well known as " missing link" in 
the AGN fueling problem. 



3.6. Impact on the detection of massive tori for SMBH growth 

In §3.2, we showed that 100 pc scale extended massive tori can be resolved and detected 
with ALMA. Here, we discuss how the detection of massive tori helps in exploring the 
SMBH growth. A plausible theoretical process of mass accretion onto a central BH is the 
turbulent viscous drag whose timescale is t v ; s ~ r 2 /z/ tu b. We adopt z/ tu b = R^rv^, where 
Rent — 100 — 1000 and are the critical Reynolds number for the onset of turbulence 
(e.g., Duschl et al. 2000; Burkert & Silk 2001) and the rotation velocity, respectively. The 
viscous time is given by t vis = -R C rit^d y n, where the dynamical timescale is determined by 
the central BH plus a surrounding massive torus system. As a result, the mass accretion 
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rate (M BH ) via the viscous drag is given as M BH oc (M t oru S /M BH ) 3/2 [l + (M BH /M torus )] 1/2 
(see also eq.(22) in Granato et al. 2004). If the massive tori with M torus /M B H ^> 1 in 
SMGs are detected with ALMA, the high mass accretion rate (close to the Eddington or 
super-Eddington accretion rate) would be a common feature of SMGs (especially hard X- 
ray detected SMGs). This statement can be checked comparing the X-ray properties (steep 
photon index, rapid variability, etc.) like those seen in NLSls in future X-ray missions. 

Moreover, KUM03 (also G04) predicts that the mass ratio M torus /M B H can be ~ 1 in 
the QSO phase, because the growth of QSO BHs has almost finished and/or we are seeing 
the final e-folding time (see Fig.l). It would be then worth comparing the massive tori in 
QSOs with those expected in SMGs with M torus /M B H 3> 1 as expected in SMGs. Thus, this 
parameter (M torus /M BH 3> 1) is a key physical conditions for the rapid growth of SMBHs 
and the appearance of QSOs and its evolution tailors the history of the BH growth process. 
Considering that the mass of a 100 pc scale dusty torus is much smaller than a central SMBH 
mass for nearby spheroidal galaxies(e.g., Djorgovski et al. 1991 and Tadhunter et al. 2003 
for Cygnus A; Jaffe et al. 1993 and Ferrarese et al. 1996 for NGC 4261; van der Marel et al. 
1993 for NGC 7052), it would be a natural interpretation that the large amount of gas around 
central BHs in SMGs accretes until exhaustion to grow to a SMBH with 10 8 — 1O 9 M on 
timescale of several 10 8 yr. Therefore, the detection of massive tori in SMGs would imply 
that the SMBH grows from ps 1O 6 M to ps 1O 9 M via the gas accretion process mainly and 
not through the merger of compact objects. 

Finally, we should comment on the case that ALMA detects less massive tori in SMGs 
than what we expect, i.e., M tovVLS / M-QYt < 1- If this is the case, it might indicate that the 
growth rate of central BHs in SMGs is much faster than the mass accretion rate from a host 
galaxy to a massive torus, or some SMGs do not have a potential to evolve into QSOs with 
M BH ~ 10 8 — 10 9 M Q . Whether or not, we suggest that exploring massive tori would cast 
light on a SMBH growth process and a formation scenario of QSOs. 

4. Discussions and Conclusions 

Given the sensitivity and angular resolution of ALMA, detailed mapping of CO emissions 
of galaxies at z ~ 2 will resolve structures on linear scales of ~ 100 pc and with mass 
-Mtoms ~ 10 8 — 10 9 M , the expected size of massive tori around SMBH. Therefore these 
capabilities can be used in order to explore the SMBH growth in SMGs, which are proto- 
spheroidal galaxies undergoing huge star formation, just before the onset of the peak of the 
QSO activity in their centre. The results of Alexander et al (2003, 2005a, & 2005b) show 
that nuclear X-ray luminosity L x > 3 x 10 43 ergs -1 is often associated to SMGs, witnessing 
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a stage in which the SMBHs have mass M BH ~ 10 6 - 1O 7 M (Borys et al. 2005). If we 
can detect massive tori with ALMA within the same population, we can definitely conclude 
that M torus /M B H ^ 1 is the key physical condition for the rapid growth of SMBHs and the 
appearance of QSOs. Moreover, the existence of massive tori in SMGs would imply that 
the SMBH grows from « 1O 6 M to ~ 1O 9 M via the gas accretion preocess rather than the 
merger of compact objects. 

In order to assess the observational feasibility, we estimate the expected number counts 
of SMGs with massive tori and check the detectability with the ALMA instrument. ALMA 
will be able to resolve and detect high-J CO emissions (J > 4) from massive tori up to 
z ~ 2 at 5(7 level in twelve hours on-source integration time and a velocity resolution of 25 
km/s. Lower abundant lines, such as HCN and HCO + lines, require much larger integration 
times at the highest angular resolution. For even high-J HCN lines at redshift 1, an on-source 
integration of at least 24 hours is necessary to detect the massive torus. Moreover, we predict 
the number count of SMGs with a massive torus (more than 10 8 M Q ) « 100 deg" 2 at the 
redshift 1 < z < 3. We stress the relevance of studying with ALMA samples of gravitationally 
lensed SMGs. Though lensing introduces the problem of lens model, nevertheless it can allow 
a much higher spatial resolution, which can be a crucial piece of information in inferring the 
gas distribution and rotation law in the very central regions around the growing SMBHs. 
Studies of lensed and unlensed SMGs would be complementary. 

It is also interesting to examine the existences of massive tori also in nearby AGNs with 
high mass accretion rate, like type I ULIRGs or NLSls. From that, it is possible to reveal 
whether the BH growth process depends on the morphology of host galaxies and redshift. 
In addition, determinations of the mass ratio of massive tori and bulges in SMGs, ULIRGs 
and NLSls would verify directly the radiation drag model. This mass ratio is basically 
determined by the energy conversion efficiency of nuclear fusion from hydrogen to helium, 
0.007 (Umemura 2001). Finally, we should emphasize that the velocity mapping on low- z 
massive tori will be crucial to reveal the formation of obscuring massive tori and also the 
fueling process from ~100 pc to ~ 1 pc. 

Acknowledgments 

We wish to thank the anonymous referee for constructive suggestions and fruitful com- 
ments. We thank M. Umemura for stimulating discussions and L. Tacconi and D. Lutz for 
helpful insights into the topic. NK appreciates that M. Umemura allows to edit his original 
figure. We acknowledge the Italian MIUR and INAF financial supports. 



-16- 



REFERENCES 

Alexander, D. M., et al. 2003, AJ, 125, 383 

Alexander, D. M., et al. 2005a, ApJ, 632, 736 

Alexander, D. M., et al. 2005b, Nature, 434, 738 

Antonucci, R., 1993, ARA&A, 31, 473 

Barvainis, R. 1987, ApJ, 412, 513 

Bahcall, J. N., et al. 1997, ApJ, 479, 642 

Bennett, C. L., et al, 2003, ApJS, 148, 1 

Boiler, Th., Brandt, W. N., & Fink, H. 1996, A&A, 305, 53 

Borys, C, et al. 2005, ApJ, 635, 853 

Burkert, A., & Silk, J., 2001, ApJ, 554, L151 

Collin, S., & Kawaguchi, T., 2004, A&A, 426, 797 

Davis, R. I., et al. 2006, ApJ, 646, 754 

Djorgovski, Weir, N., Matthews, K. & Graham, J. R., 1991, ApJ, 372, 67 
Downes, D., & Solomon, P. M., 1998, ApJ, 507, 615 
Downes D., & Solomon P. M., 2003, ApJ, 682, 37 
Dunlop, J. S., et al. 2003, MNRAS, 340, 1095 

Duschl, W. J., Stittmatter, P. A., & Biermann, P. L. 2000, A&A, 357, 1123 

Fabian, A. C, & Iwasawa K., 1999, MNRAS, 303, L34 

Ferrarese, L., & Merritt, D. 2000, ApJ, 539, L9 

Fioc, M., &Rocca-Volmerrage, B. 1997, A&A, 326, 950 

Fukue, J., Umemura, M., & Mineshige, S. 1997, PASJ, 49,673 

Ferrarese, L., Holland, F., & Walter, J. 1896, ApJ, 470, 444 

Gao, Y., et al. 2007, accepted in ApJL dastro-ph/070~3~548~l ) 



Gao, Y., & Solomon, P. M., 2004, ApJ, 606, 271 

Granato, G. L., & Danese, L., 1994, MNRAS, 268, 235 

Granato, G. L., Danese, L., & Franceschini, A. 1997, ApJ, 486, 147 

Granato, G. L., et al. 2004, ApJ, 600, 580 (G04) 

Granato, G. L., et al. 2006, MNRAS, 368, L72 

Greve, T. R., et al. 2005, MNRAS, 359, 1165 

Hao, C. N., et al. 2005, ApJ, 625, 76 

Hosokawa, T., 2002, ApJ, 576, 75 

Imanishi, M., Nakanishi, K., & Kohno, K. 2006, A J, 131, 2888 
Jaffe, W., et al. 1993, Nature, 364, 213 
Kawabe, R., et al. 1999, ASPC, 156, 45 

Kawakatu, N., Imanishi, M., & Nagao, T., 2007, ApJ, accepted dastro-ph/0702552[ ) 
Kawakatu, N., et al. 2006, ApJ, 637, 104 

Kawakatu, N., Umemura, M., & Mori M., 2003, ApJ, 583, 85 (KUM03) 
Kawakatu, N., & Umemura, M., 2002, MNRAS, 329, 572 
Kormendy, J., & Richstone, D. 1995, ARA&A, 33, 581 
Laor, A., & Draine, B. T. 1993, ApJ, 402, 441 
Laor, A., 1998, ApJ, 505, L83 
Lapi, A., et al. 2006, ApJ, 650, 42 
Magorrian, J., et al., 1998, AJ, 115, 2285 
Maiolino, R., & Rieke, G. H. 1995, ApJ, 454, 95 

Maiolino, R., Ruiz, M., Rieke, G. H., & Keller, L. D. 1995, ApJ, 446, 56 
Malkan, M. A., Gorjian, V., & Tarn, R. 1998, ApJS, 117, 25 
Marconi, A., et al. 2004, MNRAS, 351, 169 



-18- 

Marconi, A., & Hunt, L. K. 2003, 589, L21 

Mathur, S., Kurazkiewicz, J., & Czerny, B. 2000, MNRAS, 314, L17 

McLeod, K. K., & Rieke, G. H. 1995b, ApJ, 441, 96 

McLeod, K. K., & Rieke, G. H. 1995a, ApJ, 454, L17 

McLure, R. J., & Dunlop, J. S. 2002, MNRAS, 331, 795 

McLure, R. J., & Dunlop, J. S. 2001, MNRAS, 327, 199 

McLure, R. J., Dunlop, J. S., & Kukula, M. J. 2000, MNRAS, 318, 693 



Menendez-Delmestre, K., et al. 2007, ApJL, in press ( jastro-p h/0610915) 
Mineshige, S., et al. 2000, PAS J, 52, 499 

Mineshige, S., Tsuribe, T., & Umemura, M. 1998, PAS J, 50, 233 

Mochizuki, K., & Nakagawa, T. 2000, ApJ, 535, 118 

Ohsuga, K., & Umemura, M., 2001, ApJ, 559, 157 

Pounds, K. A., Done, G, & Osbore, J. P. 1995, MNRAS, 277, L5 

Radford, S. J. E., , Solomon, P. M., & Downes, D. 1991, ApJ, 369, L15 

Richstone, D., et al, 1998, Nature, 395, A14 

Salucci, P., et al. 1999, MNRAS, 307, 637 

Shankar, F., et al. 2004, MNRAS, 354, 1020 

Sato, J., Umemura, M., Sawada, K., & Matsuyama, S., 2004, MNRAS, 354, 176 
Silva, L., et al. 2005, MNRAS, 357, 1295 
Smail, I., et al. 2004, ApJ, 616, 71 

Solomon, P. M., Vanden Bout P. A, 2005, ARA&A, 43, 677 
Solomon, P. M., et al. 2003, Nature, 426, 636 



Spergel, D. N., et al, 2006, ApJ, submitted ( |astro-ph/ 0603449) 
Tacconi, L. et al. 2006, ApJ, 640, 228 



-19- 



Tadhunter, C, Marconi, A., Axon, D., Wills, K., Robinson, T. G., & Jackson, N. 2003, 
MNRAS, 342, 861 

Tielens, A. G. G. M., & Hollenbach, D. 1985, ApJ, 291, 747 Thompson, T. A., Quataert, 
E., & Murray, N. 2005, ApJ, 630, 167 

Tremaine, S., et al, 2002, ApJ, 574, 740 

Umemura, M., 2001, ApJ, 560, L29 

Umemura, M., Fukue, J., & Mineshige, S. 1997, ApJ, 479, L97 
Urry, C. M., & Padovani, P., 1995, PASP, 107, 803 



Valiante, E., et al. 2007, ApJ in press ( |astro-p h/0701816) 

Vanden Bout, P. A., 2004, ApJ, 614, L97 

van der Marel, R. P. et al. 1998, AJ, 116, 2220 

Wada, K., & Norman, C. A., 2002, ApJ, 566, L21 

Wada, K., & Tomisaka, K., 2005, ApJ, 619, 93 

Wagg, J. et al. 2005, ApJ, 634, L13 

Watabe, Y., & Umemura, M., 2005, ApJ, 618, 649 

Yu, Q., & Tremaine, S., 2002, MNRAS, 335, 965 



This preprint was prepared with the AAS IATgX macros v5.2. 



-20- 



Table. 1 The values of function H(a,/3) 
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Note. — The values in parenthesis are for 20 pc compact massive tori. The others correspond to the case of 100 pc extended 
massive tori. 

The symbol a corresponds to the fiducial case in this paper. 




Fig. 1.— A schematic sketch of the co-evolution of a SMBH and a galactic bulge (We edit 
the original figure (Fig. 1) in Umemura 2001.). The abscissa is time and the ordinate is 
arbitrary . r is the optical depth of bulges. The phase (i) denotes the early phase of a QSO 
(SMG phase). The QSO phase corresponds to the end of the phase (ii). In the phase (i), a 
massive torus around a smaller BH is predicted. 
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Fig. 2.— The required spatial resolution (arcsec) for a 20 pc scale (thin horizontal line) 
and a 100 pc scale (thick horizontal line) massive torus against redshift z. The five lines 
represent the possible spatial resolution for CO lines in the transition level (J, J — 1) with 
the ALMA instrument. 




Fig. 3.— The expected CO emission (mJy) from massive tori against redshift z in the 
fiducial case (M torus = 3 x 1O 8 M , A = 5, a = 0.5 and ft = 1). The five lines denote the 
CO flux density Sco{J, J — l)obs a t the transition level (J, J — 1) at the rest frame. The 
horizontal line represents the ALMA detection limit 0.5 mJy (50 antennas) requested 5a in 
25 km/s channel in 12 hours of on-source integration time. At the left side of filled circle, 
each line can be resolved spatially. 
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Fig. 4. — The expected HCN emission (mJy) from massive tori against redshift z in the 
fiducial case (M torus = 3 x 1O 8 M , A = 5, a = 0.5 and (3 = 1). The lines are same as Fig. 3. 




Fig. 5. — Number counts of SMGs with massive tori dn/dz [deg -2 ] against redshit z. Top 
panel shows the number counts of all SMGs. The middle one denotes dn/dz [deg" 2 ] with 
M tOTXiS > 1O 8 M . The bottom one shows the the number counts of galaxies with massive tori 
M torus > 3 x 1O 8 M . 



